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Abstract. The detection of NH3 inversion lines up to the (J, K) = (6, 6) level is reported toward the central regions 
of the nearby galaxies NGC253, Maffei2, and IC342. The observed lines are up to 406 K (for (J, K) = (6,6)) 
and 848 K (for the (9, 9) transition) above the ground state and reveal a warm (Tki n = 100 . . . 140 K) molecular 
component toward all galaxies studied. The tentatively detected (J, K) — (9, 9) line is evidence for an even 
warmer (> 400 K) component toward IC342. Toward NGC253, IC342 and Maffei2 the global beam averaged 
NH3 abundances are 1 — 2 10 -8 , while the abundance relative to warm H2 is around 10~ 7 . The temperatures and 
NH3 abundances are similar to values found for the Galactic central region. C-shocks produced in cloud-cloud 
collisions can explain kinetic temperatures and chemical abundances. In the central region of M82, however, 
the NH3 emitting gas component is comparatively cool (~ 30 K) . It must be dense (to provide sufficient NH3 
excitation) and well shielded from dissociating photons and comprises only a small fraction of the molecular 
gas mass in M 82. An important molecular component, which is warm and tenuous and characterized by a low 
ammonia abundance, can be seen mainly in CO. Photon dominated regions (PDRs) can explain both the high 
fraction of warm H2 in M82 and the observed chemical abundances. 

Key words. Galaxies: individual: NGC 253, Maffei 2, IC 342, M 82 - Galaxies: ISM - Galaxies: starburst - 
Galaxies: abundances - Radio lines: galaxies 

1. Introduction Even more difficult is the determination of local phys- 
ical parameters, such as density and kinetic temperature 
Many galaxies, including our own, show large concentra- of thfi gag Thc difficulties compared to similar studies of 
tions of molecular gas m their central few hundred parsecs. nearby douds in the galactic digk haye twQ main reasons: 
Due to the presence of bars facilitating inflow of gas, this lines of extraga i actic density or temperature trac- 
material is often originating from further out. Reaching erg; . e molecules with a high dipole moment such as CS, 
the nuclear region, the gas may trigger bursts of star for- HC3Nj Qr ^ &re brQad and weak; they require high 
mation and the formation of an active nucleus. There have quality baselines and good sen sitivity. Secondly, a low fin- 
been many studies to determine the distribution and kme- ear resolution forces us to average over regions of several 
matics of this gas, mainly using the mm-wave transitions 1QQ pc Area fining factors are nQt wdl knQwn but are cer _ 
ofCO as a tracer of the more elusive H 2 (see e.g. Combes tainly legg than unity and it becomes difficult to disentan _ 
HMH Young & Scoviflell99lJ). gle the different gas components which certainly coexist 

In the Galactic disk, there is a well established linear within one telescope beam, 
correlation between the integrated CO intensity, Iqo > and 

H2 (see Sect . 14.21 for a discussion and references). It has Nevertheless, there has been progress in determining 
become evident, however, that the relationship is different physical parameters toward the nuclear regions of nearby 
in the central regions of galaxies fe.g. Downes et al. 119921 galaxies. The most prominent examples are NGC 253, 
Sodroski et al. EHU Mauersberger et al. I1996al fl996bl M82, NGC 4945 and IC342. These and other galaxies 
Dahmen et al. I1998fl . have been mapped with high resolution in the line emis- 
sion of high density tracing molecular transitions (e.g. 
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al. 113551 Peng et al. H55Sjl . Multilevel studies in CS, HCN, 
HCO+, CH 3 CN, CH 3 CCH, CH3OH, H 2 CO and HC 3 N 
were carried out to determine H2 densities, n(H2), toward 
some of those galaxies (Mauersberger et al. 119891 11990 
1X5511 Walker et al. 115531 Jackson et al. 115551 Hiittemeister 
et al. IT5571 Paglione et al. IT5571 Seaquist et al. [15551 
Schulz et al. l2(Jl)lf) . The level populations of such molecules 
depend on their spontaneous deexcitation rates and on the 
frequency of collisions with H2, i.e. mostly on the H2 den- 
sity, and to a lesser degree on the kinetic temperature of 
the gas. From these studies of high density tracers we have 
obtained a good idea of Af(Ha) in components with n(H2) 
ranging between 10 4 and several 10 5 cm~ 3 . In order not to 
miss the less dense and more extended gas, studies of CO 
in its various transitions across the mm, sub-mm and FIR 
range were also carried out (e.g. Giisten et al. E)93s Mao 
et al. 125551 Schulz et al. !255TI Ward et al. l255Tl Bradford 
et al. 12552)1 . 

One of the less well known parameters has been the 
gas kinetic temperature. With area filling factors that are 
not known, a thermalized low density tracer like CO can- 
not be used to determine T^ n in extragalactic sources. 
Alternative thermometers for the extragalactic molecu- 
lar medium are symmetric top molecules where relative 
level populations are determined predominantly by colli- 
sions (and to a lesser degree by spontaneous radiation), 
such as NH3 or CH3CN. Temperatures in excess of 200 K 
have been obtained for clouds near our Galactic cen- 
ter (Mauersberger et al. [1986 Hiittemeister et al. [1993a 
I1993bl 1T555I Rodriguez-Fernandez et al. I255T1 Ceccarelli 
et al. l2002[l . From a few low energy ammonia lines, kinetic 
temperatures were determined toward NGC 253 (Takano 
et al. 12552)1 . Maffei2 (Henkel et al . 12555) 1. IC342 (Martin 
& HoQ]IH5J, and M 82 (Weifi et al. l2001a|> . However, sym- 
metric top molecules need H2 densities of order 10 4 cm~ 3 
to be excited. Hence they probe a dense gas component 
which is important for star formation but do not trace the 
temperature of low density, CO-emitting, molecular gas. 
For such a low density molecular component, observations 
of H2 lines are more significant. Recently, Rigopoulou et 
al. I)2()()2|l concluded from observations of rotational H2 
lines toward a number of starburst and Seyfert galaxies 
that a significant fraction (up to 10% in starburst galax- 
ies; 2 — 35% in Seyferts) of the gas in the centers of such 
galaxies has temperatures in the range of 150 K. 

In order to obtain more stringent values for the dense 
gas and to identify the origin and excitation of this gas 
we have measured higher energy inversion lines up to 
405 K above the ground state toward the central regions 
of NGC 253, Maffei2 and up to 848 K toward IC342. 
Combined with new CS J = 5 - 4 and 13 CO J = 2 — 1 
data to trace also density and column density, a dense and 
warm molecular gas component is revealed in the central 
regions of these galaxies. 



Table 1. Summary of observations 



Transition Telescope Frequency (9b c , 

GHz " 



NH 3 (1,1) 


MPIfR 100-m 


23.694 


40 


NH 3 (2,2) 


MPIfR 100-m 


23.723 


40 


NH 3 (3,3) 


MPIfR 100-m 


23.870 


40 


NH 3 (4,4) 


MPIfR 100-m 


24.139 


39 


NH 3 (5,5) 


MPIfR 100-m 


24.533 


39 


NH 3 (6,6) 


MPIfR 100-m 


25.056 


38 


NH 3 (9,9) 


MPIfR 100-m 


27.478 


34 


13 CO (2 - 1) 


HHT 10-m a 


220.399 


34 


CS (5 - 4) 


HHT 10m a 


244.936 


32 



a) in the observed frequency range, the beam efficiency of the 
HHT is 0.78 
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Fig. 1. Ammonia (NH 3 ) inversion lines (in units of mK, 
T mb ) toward the centers of NGC 253 (a 20 oo = h 47 m 33!2, 
£2000 = -25°17T6"), Maffei2 (a 20 oo = 2 h 41 m 55?2, 
S2000 = 59°36T1") and IC342 (a 20 oo = 3 h 46 m 48!6, 
£2000 = 68°05'46"), measured at Effclsberg. The spec- 
tra have been smoothed to a velocity resolution of ~ 
15 kms -1 (17 kms -1 for the (9,9) transition). Toward 
NGC 253, the feature to the left of the (1,1) line is part of 
the (2,2) profile, and the feature to the right of the (2,2) 
line is part of the (1,1) profile. 
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Table 2. NH 3 , CS(5 - 4) and 13 C0 (2-1) results 



Line 


Tmb 


J'T mb dv a 




Av a 


N(J,K) a ' b 




mK 


Kkms -1 


km s — 


km s — 


10 12 cm~ 2 


NGC253 














90 








9fURl 


NH 3 (1,1) 


41 


3.8(.8) 


288(8) 


87(22) 


50(10) 


NH,f2 2) 


19 


2 3f 7) 








NH 3 (2,2) 


29 


2.3(.6) 


303(10) 


84(18) 


23(6) 


NH 3 (3,3) 


65 


5.1(.6) 


177(3) 


72(6) 


45(5) 


NH 3 (3,3) 


63 


9.7(.8) 


295(5) 


145(13) 


84(7) 


NH 3 (4,4) C 


14 


1.4(.4) 


325(13) 


95(23) 


11.2(3.1) 


NH 3 (6,6) 


20 


1.7(.5) 


169 


79(20) 


12.9(3.8) 


NH 3 (6,6) 


21 


1.2(.4) 


303 


55(18) 


9.2(3.0) 


CS(5-4) 


45 


7.7(0.9) 


215(10) 


162(20) 




13 CO(2- 1) 


244 


52.4(4.3) 


251(8.5) 


202(19) 




Maffei 2 












NH 3 (6,6) 


3.4 


.21(.06) 


-74(8) 


59(17) 


1.5 (.4) 


NH 3 (6,6) 


5.1 


.39(.06) 


10(6) 


72(12) 


2.8 (.4) 


CS(5-4) 




< 0.5 d 








13 CO(2- 1) 


6.2 


5.2(1.6) 


-99(10) 


79(16) 




13 CO(2- 1) 


6.2 


8.7(1.9) 


23(14) 


132(32) 




IC 342 












NHs(l,l) 


26.3 


1.5(.07) 


31.4(1.2) 


52.8(2.8) 


20.0(1.0) 


NH 3 (2,2) 


19.4 


1.7(.10) 


26.8(2.2) 


83.0(6.7) 


16.7(1.0) 


NH 3 (3,3) 


26.4 


1.7(.06) 


35.1(1.1) 


58.8(2.3) 


14.9(.5) 


NH 3 (4,4) 


8.8 


.47(.04) 


29.1(2.2) 


48.9(4.7) 


3.8(.3) 


NH 3 (5,5) 


5.6 


.29(.03) 


34.8(2.6) 


47.6(6.2) 


2.2(.2) 


NH 3 (6,6) 


7.0 


.45(.06) 


37.4(3.7) 


60.7(8) 


3.2(0.4) 


NH 3 (9,9) 


5.8 


.34(.08) 


35 


60 


2.1(.5) 


CS(5 - 4) 


8.5 


.78(.17) 


30 


86(21) 




13 CO(2- 1) 


253 


17.2(1.1) 


24.5(1.9) 


64.1(4.5) 




M82 












CS(5-4) 


9 


1.3(0.3) 


150(11) 


129(36) 




CS(5-4) 


17 


1.8(0.2) 


303(6) 


94(13) 




13 CO(2- 1) 


100 


13.5(1.2) 


150 


129 




13 CO(2- 1) 


150 


15.9(1.0) 


303 


94 





a) values in italics have been fixed in the Gaussian fits, all errors are lcr 

b) total beam averaged column density of the observed (J,K) rotational level including both the upper and lower inversion 
state. 

c) (4,4) data toward NGC 253 may be affected by pointing problems (see Sect. l3.ll 

d) lcr upper limit 



2. Observations 

2.1. EfFelsberg Observations of NH3 lines 

Between October 2000 and September 2002, the (J, K) = 
(1,1) to (6,6) inversion lines of ammonia were observed 
using the EfFelsberg 100-m telescope 1 . The observed line 
frequencies and the corresponding beam sizes are summa- 
rized in Tabled 



1 The 100-m telescope at Effelsberg is operated by the Max- 
Planck-lnstitut fur Radioastronomie (MPIfR) on behalf of the 
Max-Planck-Gesellschaft 



The data were recorded with a dual channel K-band 
HEMT receiver with T sys ~ 200 K on a main beam bright- 
ness temperature scale (T m b). The measurements were 
carried out in a dual beam switching mode, with a switch- 
ing frequency of 1 Hz and a beam throw of 2' in azimuth. 
An 'AK90' autocorrelator included eight spectrometers 
with 512 or 256 channels and bandwidths of 40 or 80 MHz 
each, leading to channel spacings of ^1 or 4kms _1 . The 
eight spectrometers were configured such that the (J, K) 
= (1,1) to (4,4), (4,4) and (5,5), or (5,5) and (6,6) lines 
could be observed simultaneously. Data with one individ- 
ual line, either (J,K) = (5,5) or (6,6), were also taken. 
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Fig. 2. CS J=5-4 and 13 CO J=2-l profiles measured 
with the HHT. For a comparison with NH3 data from 
Wcifi et al. f2001a)l we also include profiles measured to- 
ward M82 (a 20 oo = 9 h 55 m 52!6, <5 2 ooo = 69°40'46"). The 
velocity resolution has been smoothed to 17 kms -1 . 



At least two backends were centered on each ammonia 
line to sample both linear polarizations independently. For 
NGC 253, only one of the two receiver channels provided 
useful data. 

The (9,9) line at 27.478 GHz was observed in December 
1995 using a single channel 1 cm HEMT receiver equipped 
with a 1024 channel autocorrelator. System temperatures 
were near 500 K on a T m b scale. The beam size was 34". 
The (9,9) line was measured in a position switching mode 
with a bandwidth of 25 MHz and a channel spacing of 
-0.27 km s -1 . 

Flux calibration was obtained by regularly observing 
both the continuum and the ammonia lines of W3(OH) 
(for fluxes see Mauersberger et al. 119881 Ott et al. H994|) . 
Pointing was checked every hour. W3(OH) was the point- 
ing source for IC 342 and Maffei 2; PKS 0023 - 26 was used 
for NGC 253. The pointing accuracy was found to be sta- 
ble within 5"-10". 

A linear baseline was removed from each individual 
spectrum and spectra were added with relative weights 



depending on the noise level outside the velocity range of 
the respective line. We estimate the flux calibration of the 
final reduced spectra to be accurate within ±15%. 

2.2. CS and 13 CO observations with the 
Heinrich-Hertz- Telescope 

The J=5-4 transition of 12 C 32 S (hereafter CS) and 
the J=2-l transition of 13 C 16 (hereafter 13 CO) 
were observed in November 2001 toward NGC 253, 
Maffei 2, IC342, and M 82 using the 10-m Heinrich-Hertz- 
Telescope 2 (Baars & Martin H996[) . At the line frequencies 
of 244.9536GHz (CS) and 220.3987GHz ( 13 CO), the beam 
width was —32". The zenith optical depth at 225 GHz as 
determined by a tipping radiometer was typically 0.2 — 
0.4. A single channel SIS receiver with a total bandwidth 
of —0.8 GHz tuned in double sideband mode was used for 
the observations. The spectrometer was a 2048 channel 
AOS with a channel separation of 0.48 MHz (—0.6 kins" 1 ). 
The beam efficiency of the HHT is 0.78 in the observed 
frequency range. System temperatures were 600 — 1000 K 
on a T m b scale. The calibration procedure was described 
by Mauersberger et al. I|1999|) . The observing mode was 
symmetric wobbler switching with a throw of ±200" in 
azimuth. The pointing was checked every 3 — 5 hours and 
was accurate to <5". Errors in the calibration procedure, 
atmospheric fluctuations and deviations from a sideband 
gain ratio of unity should result in an uncertainty of ±20% 
in the intensity scale. After co-adding the spectra, we sub- 
tracted first order baselines and smoothed to a velocity 
resolution of — 17 kms -1 . 



3. Results 

3.1. NH3 

The measured spectra are displayed in Fig.^ and the re- 
sults of Gaussian fits to the profiles, namely the velocity 
integrated main-beam brightness temperatures J T m bdv, 
the velocities relative to the local standard of rest (i>lsr) 
and the full- width to half power linewidths (Avi/2) are 
summarized in Tabled For the (9,9) line toward IC 342 
we have fixed the central velocity and lincwidth in order 
to reduce the nominal error for the integrated intensity. 
We have clearly detected the (1,1) . . . (6,6) transitions 
and tentatively also the (9,9) transition toward IC342, 
the (1,1) . . . (4,4) and (6,6) lines toward NGC 253 and the 
(6,6) line toward Maffei 2. These are the first detection 
of the (5,5) transition, the first definite detections of the 
(6,6) transition and the first tentative detection of an NH3 
(9,9) line outside the Galaxy. 

Toward NGC 253, we find two velocity components 
in the (1,1), (2,2), (3,3) and (6,6) lines at -175 and 
290 kms -1 , each with a linewidth of 80-100 kms -1 , which 
presumably arise from two intensity peaks also seen 



2 The HHT is operated by the Submillimeter Telescope 
Observatory on behalf of Steward Observatory and the MPIfR. 
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in higher resolution maps (Mauersberger et al. 2996a 
Harrison et al. 1999). These velocity components cannot 
be distinguished in our (4,4) profile. As line width and cen- 
tral velocity indicate, this line only shows the 300 kms -1 
component. Therefore we cannot exclude an exception- 
ally large pointing error favoring the higher velocity peak. 
From CO mapping (e.g. Mauersberger et al. 1996a, Dumke 
et al. I2001f) the observed line profile can be reproduced if 
the telescope was mispointed by ~ 10" toward the SW. 
Therefore we will not account for the detailed shape of 
the (4,4) line in the further interpretation of our results 
but will take the presence of the (4,4) line as just another 
indication for the presence of warm molecular gas. 

For the lower lying transitions of ammonia and other 
molecules toward Maffei2, there are two velocity compo- 
nents, at about -80 and +6 km (e.g. Henkel et al. EOOO ), 
each with a line width of ~ 50 kms -1 . Our (6,6) spectrum 
is consistent with these results (see Fig.^and Table 1). A 
single velocity component could fit all data from IC342. 



3.2. CS(5-4) and 13 CO (2-1) 

We have detected CS (5 - 4) emission toward NGC253, 
IC342, and M82. Toward Maffei2, we obtain an upper 
limit. 13 CO emission was observed toward all four galaxies. 
The profiles are shown in Fig.Eland results from Gaussian 
fits are presented in TableEl 



Table 3. CS peak line intensities in main beam brightness 
temperature (mK) at 32" resolution, the HHT CS J =5-4 
beam size. Estimated errors are ±30% (see Sect. 13. 3j) . 



Transition NGC 253 IC342 M82 (SW) 



CS J = 


2 




1 


160" 


80 a 


92 d 


CS J = 


3 




2 


77° 


45 a, e 


34 a 


CS J = 


5 




4 


50 a ,45 i ' 


12 tt , 8.5 6 


l3 b,e 


C 34 S J 




2 


- 1 


23 a 


12 c ' f 


14« 


C 34 S J 




3 


- 2 




8 c,f 


4 a 



a) original data from Mauersberger & Henkel ( 1989 ) 

b) this paper 



c) original data from Mauersberger et al. ( 1995 I 

d) original data from Baan et al. ( 1989 1 

e) data measured with 0", +10" offset to other lines 

f) data measured with 0", —10" offset to other lines 



4. Discussion 

In the following, we discuss the nuclear regions of all 
four galaxies observed in the 13 CO J=2— 1 and CS J=5-4 
transitions, i.e. NGC 253, Maffei 2, IC342, and M82 (see 
Fig.EJ). NH 3 data from M 82 will be taken from Weifi et al. 
( 2001a). The NH 3 observations (Fig.Q] Weifi et al. l2001afl 
will be used to derive kinetic temperatures in regions dense 
enough (i> 10 4 cm -3 ) to excite the NH3 inversion lines. 
The 13 CO observations are needed to estimate H2 column 
densities and the CS data will be used to determine n(H 2 ) 
densities, CS column densities, and CS abundances rela- 
tive to H2. 



3.3. Source structure, beam size and line intensities 

To constrain molecular abundances and physical condi- 
tions from observed CS line intensities, it is crucial that 
the transitions are observed with comparable beam sizes 
or that appropriate corrections are applied accounting for 
beam size and source morphology. Not included in this 
discussion is Maffei 2 for which not enough CS data exist. 
For IC 342 we have assumed that the observed CS and 
NH3 emission arises from the same volume as the 13 CO 
J=l— emission that was observed with high spatial reso- 
lution by Meier et al. (HEP) . Smoothing the 13 CO J=l-0 
data cube to different beam sizes thus enables us to in- 
vestigate the dependence of observed line intensity on the 
beam size. The same procedure was applied to the distri- 
bution of the molecular gas in M 82 using the high angular 
resolution C 18 J=l-0 map of Weiss et al. (|2001bf) and to 
NGC 253, using the distribution of CS J=2-l (Peng et al. 
1996). The CS line intensities expected at 32" resolution 
(the beam size of our CS J=5-4 data, see Sect. 2.2) are 
summarized in TableOl Given the uncertainties of the lat- 
ter approximation and accounting for possible differences 
between the C 18 or 13 CO versus the CS distributions in 
M82 and IC342, we estimate an error of ±30% for the 
intensity of individual CS lines convolved to a beam size 
of 32". 



4.1. A//-/3 column densities and rotational temperatures 

In the optically thin case, the beam averaged column den- 
sity of a rotational state of ammonia (i.e. upper plus lower 
inversion level) is related to the velocity integrated main- 
beam brightness temperature by 

1.5510 14 cm- 2 J(J+l) f 
N(J,K) = — — / T mb dw, (1) 

with v being the line frequency in GHz, and J T m \,dv 
being the integrated intensity in Kkms -1 . Here we as- 
sumed that the excitation temperature of the inversion 
lines is > IK, the energy difference of the two states of a 
given (J, K) level. The calculated total column densities 
and their errors are displayed in Table Because most 
lines were measured simultaneously, the main contribu- 
tion to the given uncertainties arises from the errors in 
the Gaussian fits. 

If one plots \og 1Q (N(J,K)/(g op (2J+l))), where g op = 
2 for ortho-ammonia [K = 3,6,9,...) and 1 for para- 
ammonia (K = 1,2,4,5,7...), against the energy of 
the involved (J, K) rotational state, the slope of a lin- 
ear fit to two or more rotational levels, m, is related 
to a rotational temperature (T rot ) by T rot = — 0.434/m 
(e.g. Mauersberger et al. H986|l . Such rotational diagrams 
(Boltzmann plots) are displayed in Fig.|2|for the observed 
galaxies. For NGC 253, a Boltzmann plot is given for each 
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Fig. 3. A Boltzmann plot for the metastable NH3 levels (i.e. with J = K) observed toward Maffei 2 (data for (J, K) < 
(6,6) are from Henkel et al. [2000 ) , for IC342 and for the two velocity components of NGC253 at 170 and 300 
kms -1 . We also show the M82 data from Weifi et al. ( 20 01afl correcting an error in the position of the (4,4) limit in 
their paper. The x-axis denotes the energy level of the corresponding rotational transition in K, the y-axis denotes 
log(N(J,K))/g op (2J+l) (see Sect.EJ. For abetter comparison of the resulting temperatures, the scales in all frames 
and axes are the same. Linear regressions to the data and the resulting values for T rot are shown. 



of the two observed velocity components. Column den- 
sities for Maffei 2 have been calculated from the data in 
Henkel et al. <|2UUUjl . 

Fig-El shows the first extragalactic Boltzmann plots 
from ammonia with fits extending beyond the (4,4) line 
(energy level: 199 K) to energies of 405 K and 848 K. 
Tabled summarizes data for the central regions of the 
galaxies discussed in more detail below. Model calcula- 
tions (e.g. Walmsley & Ungerechts 1983) indicate that 
due to collisional deexcitation of metastable levels (i.e. 
with J = K) via lower lying metastable levels the rota- 
tional temperature is only a lower limit to the gas kinetic 
temperature. This effect is less pronounced for the higher 
metastable ammonia lines, which therefore should be a 
more reliable indicator of kinetic temperatures. 

4.1.1. Comments on individual sources 

NGC253. For NGC253, the (4,4) line has not been in- 
cluded in the fit (see Sect. l3.lf) . We show fits for two 
velocity components, each corresponding to one of the 
twin peaks separated symmetrically by about 10" from 
the center of NGC253 (Mauersberger et al. H996afl . The 
data can be fitted by a single rotational temperature of 
142(±14)K for the 170 kms^ 1 component and 100(±3)K 
for the 300 km s -1 component. These values are much 



higher than the rotational temperature of 17-50 K derived 
from measurements of the (1, 1) to (3, 3) lines and a non- 
detection of the (4,4) line by Takano et al. $M^ . Ortho 
to para ammonia abundance ratios of up to 6 have been 
claimed by Takano et al. 1)2002(1 for some velocity compo- 
nents toward NGC253. Such a high ratio would indicate 
that NH3 was formed in a very cold (< 10 K) environment. 
Although the (3,3) line seems to be stronger than expected 
from the intensities of the (2,2) and (4,4) lines, the (6,6) 
line (which also belongs to the ortho species) has a nor- 
mal intensity. We therefore cannot confirm an anomalous 
ortho-to-para ratio in NGC253. Unlike in Maffei 2 and 
IC 342 (see below), in NGC 253 we cannot identify an im- 
portant amount of cool (< 50 K) in ammonia. In Sect l4.4l 
we argue that this is due to an elevated relative abundance 
of NH3 in the warm component. 



If we assume a thermalization at those temperatures 
and further assume that the non-metastable ammonia lev- 
els (i.e. those with J > K) are not populated (a reasonable 
assumption because a significant population would require 
n(H2) > 10 6 cm~ 3 , e.g. Mauersberger et al. H985f) . the to- 
tal ammonia column densities obtained are 1.610 14 cm -2 
for the 170 kms- 1 and 2.110 14 cm~ 2 for the 300 kms" 1 
component. 
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Maffei2. As toward NGC253, two velocity components 
are detected in Maffei 2. While Henkel et al. I|2000[l sug- 
gested that the -80kms _1 component shows higher exci- 
tation than the +6kms~ 1 feature, our (6,6) line spectrum 
(Fig.^) indicates equal excitation within the limits of ob- 
servational accuracy. A single rotational temperature does 
not yield a good fit to the data. The (1,1) and (2,2) lines 
can be fitted by T rot = 48(±15) K and the (3,3) to (6,6) 
lines by T lot = 132(±12)K. In a previous study, Takano 
et al. I|2()()()J) derived from a stronger than expected (3,3) 
line of ortho- ammonia an anomalously high ortho/para 
ratio of 2.6 for ammonia, which would be an indication of 
NH3 formation in a cold (13 K) environment. Our higher 
signal-to-noise (4,4) and (6,6) data do not support this 
finding: the ortho/para ratio toward Maffei 2 is normal. 
Assuming that the (0,0) rotational level is connected to 
the (1,1) and (2,2) levels by a rotational temperature of 
48 K, the total ammonia column density in all metastable 
levels is 8.6(1.3) 10 13 cm" 2 . 

IC342. There is no single rotational temperature for all 
the observed lines. The (1,1)... (3, 3) transitions can be 
represented by T rot = 53(±1) K and the (5, 5) . . . (9, 9) 
lines by T rot = 443(±130) K. If we had not included the 
(9,9) data, the temperatures obtained would have been 
similar to those obtained for NGC 253. Assuming that the 
(0,0) rotational level is connected to the (1, 1) . . . (3, 3) 
levels by a rotational temperature of 53 K, the to- 
tal ammonia column density in all metastable levels is 
7.5(0.5) 10 13 cm -2 . About 1/3 of this total column belongs 
to a hot component of several 100 K. 

M 82. Although we are adding no new ammonia data to 
those observed by Weifi et al. (2001a), we notice that in 
their Boltzmann plot the limit for the (4,4) line was er- 
roneously plotted too low. While the (1,1) . . . (3,3) lines 
can be described by a common temperature of 29(±1) K, 
the (4,4) line's upper limit cannot exclude that the higher 
ammonia levels are described by components with high 
rotational temperature, like in other galaxies observed. 

4.2. Determining N(R 2 ) from 13 CO 

The integrated intensity of 12 CO 1 — has been empiri- 
cally related to the H2 column density for Galactic disk 
molecular clouds (e.g. Dickman 119751 Sanders et al. 119851 
Blocmen et al. I1986[l . The theoretical explanation of such 
a universal X-factor requires three conditions to be ful- 
filled (Dickman et al.[1986J Mauersberger & Henkel 1993), 
namely that (a) there are many unresolved clouds with 
optically thick CO emission within a beam (CO counts 
clouds), (b) that these clouds are in virial equilibrium 
and (c) that the CO rotational levels are populated via 
collisions with H 2 molecules. For molecular clouds near 
galactic centers (including our own), conditions a) and 
c) are likely fulfilled. It is, however, by no means clear 
whether molecular clouds close to the central regions of 



galaxies are virialized (except perhaps the very densest re- 
gions). Tidal forces from the nuclear regions or high stellar 
densities may easily rival the proper gravitational forces 
of a molecular cloud near a galactic center (e.g. Glisten 
1989, Mauersberger et al. I1996a|l . For such tidally dis- 
turbed clouds, the theoretical explanation of the X factor, 
which works so well for molecular clouds in the Galactic 
disk, would not be valid anymore. This is even more true 
for any molecular intercloud medium (e.g. Downes et al. 
nMSl Jog & Das ITMfl Solomon et al. ITM71 Downes & 
Solomon 1998), which is certainly not virialized. An in- 
creasing number of examples reveal that the same X factor 
from the Galactic disk cannot be used near the centers of 
galaxies (Galactic Center: Sodroski et al. 1 19951 Dahmen et 
aLEm NGC253: Mauers berger et a,l. H996al NG C 4945: 
M auersbe rger et al. H996bl M82: Mao et al. 12113 Weifi et 
al. l2001bl ultraluminous galaxies: Solomon et al. [1997 ) . 

Relating optically thin transitions of CO to the H2 
column density requires far fewer assumptions than for 
optically thick 12 CO transitions. The mm-transitions of 
12 C 18 (hereafter C 18 0) are in most conceivable cases 
of astrophysical interest optically thin (e.g. Mauersberger 
et al. I1992|l : the transitions of 13 CO should be in most 
cases optically thin or only moderately optically thick. 
Because of the high chemical stability of CO and its iso- 
topes their relative chemical abundance should vary less 
than that of other molecules in the interstellar medium. 
Model calculations (Mauersberger et al. I1992|) show that 
for a large range of temperatures and densities the in- 
tensity / = J T m bdu (in Kkms -1 ) of an optically thin 
isotopic CO ('C^O) J = 2 - 1 line is related to the H 2 
column density by 

iV(H 2 ) = 5.3 10 14 cm- 2 (^^)- 1 /( i C J O). (2) 

[H 2 J 

Any increase in T cx of the corresponding levels due to 
higher excitation (i.e. higher n(H 2 ), 7ki n ) is approximately 
compensated by a loss of level population of the corre- 
sponding states to higher states in the J ladder; this is why 
Eq.[21holds even if the J — (2 — 1) line is subthermally ex- 
cited. An unknown quantity in Eq.|21 is the relative abun- 
dance of the CO isotopomers studied. For Galactic disk 
molecular clouds, a CO abundance of 8 10 -5 is frequently 
assumed (Frerking et al. I1982D . Toward the centers of spi- 
ral galaxies, metallicities are typically higher by a factor 
of two than in the local disk (see Vila-Costas & Edmunds 
1992}. On the other hand in the centers of galaxies the 
amounts of carbon in the form of C I and CO may be 
comparable (Gerin & Phillips I2UUU1 Israel & Baas I2UU2"|) . 
It is therefore plausible to expect for the galaxies studied 
here a global CO/H 2 abundance similar to the Galactic 
disk value, with an uncertainty which could be a factor of 
two in either direction. 

Model calculations suggest that CO isotopomers are 
not very greatly affected by selective dissociation (Chu & 
Watson 1983, Glassgold et al. I1985|) . A fractionation of 
13 CO via the reaction 

12 CO + 13 CO+ -^ 13 CO + 12 C + AE 35K (3) 
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(Watson I1976() becomes significant only for kinetic tem- 
peratures smaller than 35 K (e.g. Langer et al.[l984). As 
long as the kinetic temperatures of the CO emitting gas 
are higher than this it is safe to assume that fractiona- 
tion of 13 CO introduces an error of the H2 column density 
which is negligible in comparison to the uncertainty of the 
CO/H2 abundance. 

The 12 C/ 13 C and 16 0/ 18 ratios have been deter- 
mined to b e 50 an d 200 (Henkel & Mauersberger 119931 
Weifi et al. 1200 11 it in the galaxies studied. Eq. (1) trans- 
forms into AT(H 2 ) = 1.3 10 21 cm" 2 7(C 18 0) and 7V(H 2 ) = 
3.3 1 20 ciTr 2 /( 13 CO). We have used the latter equation 
together with our measurements of the 13 CO ,7 = 2—1 
line to determine the beam averaged H2 column densi- 
ties summarized in Table0] If the 13 CO lines are slightly 
optically thick, the iV(H 2 ) values may be higher (see be- 
low). While C 18 lines could be an even more accurate 
tracer of ./V(H 2 ), such lines are substantially weaker and 
are difficult to detect with a proper signal-to-noise ratio. 

A comparison between the derived H2 column den- 
sities and 12 CO J = 1 — measurements (account- 
ing for different spatial resolutions, see Sect l3.3l) to- 
wards NGC253 (Houghton et al. HMTl Sorai et al. l2TTO 
and IC342 (Crosthwaite et al. IMiTl Meier & Turner 
2001) shows, that a conversion factor Xqo which is 
a factor of 3 — 5 smaller than the Galactic value of 
2.3 10 20 cur 2 (Kkms -1 )- 1 (Strong et al. lT3gg|> is appro- 
priate for starburst nuclei. This is in line with recent de- 
terminations of Xqo in M82 (Weifi et al. l2001b|> . 

4.3. n(R2) and CS abundance 

We now use the intensity of the CS (5 — 4) lines 
given in Tabled to complement the multilevel studies in 
Mauersberger & Henkel l|1989J) . The levels corresponding 
to this transition need an H 2 density ~ 10 4 cm~ 3 to be 
excited. We used the large velocity gradient (LVG) code 
described in Mauersberger & Henkel fiJBgt to investigate 
the H 2 density, the CS column density and relative abun- 
dance (n(H 2 ), AT(CS) and A"(CS)=JV(CS)/JV(H 2 )). Note 
that the values given below for AT(CS) are beam aver- 
aged values for the CS emitting dense gas component (the 
LVG code gives the source averaged iV(CS) which can be 
multiplied with the beam filling factor, given by a com- 
parison of the code's results with the observed intensi- 
ties, to obtain the beam averaged iV(CS)). Any lack of 
CS emission in the intercloud component due to a low 
abundance or insufficient excitation increases the local 
abundance in the CS emitting component. LVG models 
were calculated for log(n(H 2 )) = 2.0 to 6.0 in steps of 
0.2 and log(AT(CS)/grad(i;)) = -11.7 to -6.7 in steps of 
0.2. The relative abundance [CS]/[C 34 S] was fixed to 25 , 
the Galactic 32 S/ 34 S abundance ratio (Chin et al. lTM)|) . 
Since the kinetic temperature of the gas is high (see the 
results in Sect. 14.1(1 and the excitation of CS does only 
weakly depend on T kin for temperatures above 30 K (see 
e.g. Mauersberger & Henkel we fixed T kin to 60 K. 



Solutions were selected by comparing all beam corrected 
CS and C 34 S line intensity ratios derived from Tablo|21to 
the corresponding LVG line intensity ratios using a % 2 - 
test. Only those solutions were permitted where the beam 
averaged H 2 column densities were not below the values 
derived from our 13 CO J — 2 — 1 measurements (NGC 253 
and IC342). 

For Maffei2, our non-detection of the J=5-4 line 
and the intensity of the J=3-2 CS line detected by 
Mauersberger et al. (1989a) do not permit a reliable de- 
termination of n(H2). For IC 342 all observed line intensi- 
ties can be fitted simultaneously by the LVG models. We 
derive a moderate n(H2) density of 10 3 ' 9±0 ' 3 cm -3 . The 
beam averaged CS column density is between 1 10 13 cm~ 2 
and 610 13 cm -2 . Using the H2 column density given in 
Tabled this corresponds to a beam averaged relative CS 
abundance of i[p 8 - 4±0 - 5 . From an HCN multilevel study 
Paglione et al 1(1997(1 estimated elevated kinetic temper- 
atures and densities toward IC342 which are consistent 
with our models. Using more transitions of HCN, Schulz 
et al. 1)2001(1 point out that 5% of the dense gas might 
have densities of 10 6 cm -3 and might be as cool as 30 K. 

For NGC 253 and M 82 not all observed line intensi- 
ties can be fitted by the LVG models simultaneously. As 
discussed in Mauersberger & Henkel lfK?55jl for NGC 253 
this indicates, that the lower transitions of CS arise from a 
lower density component containing dense concentrations 
traced by the 5 — 4 transition. This explanation is sup- 
ported by the LVG solutions we obtain when we do not 
include either the CS J = 2 — 1 or CS J = 3 — 2 transition: 
for NGC 253 we obtain n(H 2 ) = 10 4 2±0 - 7 cm" 3 (using CS 
J = 2 - 1, CS J = 5 - 4 and C 34 S J = 2 — 1) and 
n(H 2 ) = 10 4 - 5±a8 cm- 3 (using CSJ = 3-2, CSJ = 5-4 
and C 34 S J = 2 — 1). Including only the low J transitions 
yields n(H 2 )) = 10 3 7±0 - 6 cm -3 . The beam averaged CS 
column density is between 210 13 cm~ 2 and 210 14 cm -2 . 
The beam averaged CS abundance therefore is io -8 - 4±0 - 5 . 
Our density estimates are consistent with multilevel stud- 
ies of HCN (Paglione et al. ITWjl . HC 3 N (Mauersberger 
et aL HTTO . CH 3 CN and CH 3 CCH (Mauersberger et al. 

nmn . 

The corresponding results for M82 are: n(H 2 )) = 
10 3.9±o.7 cm -3 ( using CS J = 2 — 1, CS J = 5 — 4, C 34 S 
J = 2 - 1, and C 34 S J = 3 - 2) and n(H 2 )) = 10 4 2±0 - 3 
(using CS J = 3 - 2, CS J = 5 - 4 and C 34 S, J = 2 - 1 
and C 34 S J = 3 — 2). The density is not well constrained 
using the low J lines only. LVG solutions then cover a den- 
sity range from n(H 2 ) = 10 3 3 to 10 5 5 cm -3 and low H2 
densities are associated with high 'local' CS abundances 
and vice versa. Interestingly, solutions are only permit- 
ted for Tkin < 25 K. Beam averaged CS column densities 
including the CS J = 5 — 4 line are in the range be- 
tween 110 13 cm -2 and 610 14 cm~ 2 . The beam averaged 
CS abundance is i[p 7 - 9±0 - 8 . The LVG results are sum- 
marized in Tabled Our density estimates are consistent 
with those from HCN (Paglione et al. lT59T| and CH 3 CCH 
(Mauersberger et al. 11991(1 . 
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Table 4. Physical parameters in the central regions (~ 35") of nearby galaxies 



Source 


D 


35" equiv. 


iV(H 2 ) a 


logn(H 2 ) 


Trot 


iV(NH 3 ) 


iV(CS) 


logpf(NH 3 )) 


log(X(CS)) 






beam size 


















Mpc 


pc 


10 21 cm" 2 


cm -3 


K 


10 13 


cm -2 






NGC 253 


2.5" 


420 


17.4 


4.2 ± 1.1 


100. . . 142 


37(7) 


2-20 


-7.7(0.3) 


-8.4 ±0.5 


Maffei 2 


2.5 C 


420 


4.6 




48, 132 


8.6(1.3) 




-7.7(0.3) 




IC 342 


1.8 d 


300 


5.7 


3.9 ±0.3 


53, 443 


7.5(.5) 


1-6 


-7.9(0.3) 


-8.4 ±0.5 


M82 


3.9 e 


660 


9.7 


4.3 ± 1.0 


29 ; 


1.0(.3) / 


1-60 


-9.0(0.3) 


-7.9 ±0.8 



a) from our 13 CO J = 2 — 1 data; the uncertainty is a factor of 2 in both directions, b) Mauersberger et al. H996al c) Karachentsev 
et al. rtWl d) McCallEjlH e) Sakai & Madore H*9"9l?l f) Weifi et al. I2001bl 



Table 5. Masses and NH 3 abundances of the warm (> 
150 K) gas component 





M a 

iU warm 

10 6 


A^total 

M 




logX(NH 3 ) c 


NGC 253 


2.5 


75 


0.03 


-6.1 


IC342 


1.3 


13 


0.10 


-7.5 


M82 


75 d 


100 


0.2 . . .0.8 d 


< -9 


GC e 






0.3 


£-7 



a) Within a 14" x 27" aperture, from Rigopoulou et al ( 2002 ) 
scaled to the distances in Tabled 

b) from 13 CO data in this paper, the uncertainty is a factor 



of 2 

c) The warm NH3 relative to the warm H2 

d) It is unknown to us to which pointing position and 
orientation of the ISO beam in Rigopoulou et al. ( 2002 ) the 
mass in Col. 2 refers; therefore the large uncertainty of Col. 4 

e) Rodriguez-Fernandez et al. 12001H 



4.4. Origin and heating of the warm gas 

In the complex environment of the central few 100 pc 
of an active galaxy there are certainly several coexist- 
ing gas components. Transitions from molecules with a 
high dipole moment such as ammonia and CS need a den- 
sity ~ 10 4 cm -3 to be excited, in contrast to CO lines 
which need a few 100cm -3 . Our ammonia Boltzmann 
plots also suggest that there are several components with 
different temperatures within our beams. This is also 
what we observe toward the central region of our Galaxy 
(Hiittemeister et al. I1993b(l . The relat ive abundances of 
the high density tracers NH3 and CS in Tablc^] are from a 
comparison with the total molecular column densities de- 
rived from 13 CO; they are therefore an average over all gas 
components. We can, however, estimate the relative abun- 
dance of warm ammonia in the gas component where it 
probably arises. This component is traced by rotational 
quadrupole transitions of H 2 , e.g. the — S(0) line, 
the upper level of which is 510 K above the ground state. 
From such lines masses of warm gas with a temperature of 
~ 150 K have been recently determined by Rigopoulou et 
al. (2002) for a number of starburst and Seyfert galaxies. 
We have scaled the masses quoted in Rigopoulou et al. 
(2005J) to the distances used here and list them in Table 
|5J When we compare the masses of warm H2 toward M 82, 
IC 342, NGC 253 with the total H 2 masses from our 13 CO 



data and also take into account Galactic Center H2 mea- 
surements obtained by Rodriguez- Fernandez et al. (2001) 
it turns out that the warm gas makes up a substantial 
fraction of the total gas mass in the centers of nearby 
starburst galaxies and our own Galaxy. We cannot con- 
firm, however, the claim by Bradford et al. |(2D02) that 
all the gas traced by CO is warm (and, hence, requires a 
uniformly heating source such as cosmic rays) since this 
would lead to contradicting H2 mass estimates from ob- 
servations of H2 and CO. 

Newly formed high mass stars may heat only very lo- 
calized regions via dust absorption. Even in the central 
regions of luminous galaxies, global dust temperatures 
rarely exceed 50 K (e.g. Henkel et al. 119861 1. In the case of 
NGC 253, the central mid-infrared emission comes mainly 
from cold (23 K) dust; only a small fraction, <C 20% of 
the dust mass has temperatures similar to 148 K (Melo 
et al. 12002(1 . Possible mechanisms capable of elevating the 
H2 temperature to global values of 100 K and more over 
an extent of several 100 pc are shocks, which may result 
e.g. from the dissipation of tidal motions, photodissoci- 
ation and irradiation by X-rays (see e.g. Rigopoulou et 
al. 1200211 . Other possible heating mechanisms involve cos- 
mic rays (Glisten et al. 119811 Farquhar et al. 1 1994(1 or a 
magnetic ion-slip process (Scalo 1977). 

In the central regions of IC 342 and NGC 253, the rel- 
ative abundance of the hot ammonia to the hot H2 com- 
ponent is elevated over the relative ammonia abundances 
observed also including the cold component. In NGC 253 a 
value of 10~ 6 is reached, which is similar to values found in 
Galactic hot cores (e.g. Hiittemeister et al. H993b| and the 
central region of the Milky Way (Rodriguez-Fernandez et 
al. 12001(1 . The high NH 3 abundances in these galaxies and 
presumably also in Maffei 2 exclude heating mechanisms, 
such as photo-ejection from grains (Genzel et al. [T984), 
capable of destroying this fragile molecule which has a 
dissociation energy of only 4eV. NH3 has been observed, 
with a similar abundance as in the warm components of 
NGC 253 and IC342, in Wolf Rayet environments (Rizzo 
et al. I2001() . which can be explained by a release of NH3 
from dust grains facilitated by weak shocks from the ex- 
pansion of the bubble created by the stellar wind. C-shocks 
have also to be used to explain the abundance of com- 
plex molecules and the heating in the central region of our 
Galaxy (Flower et al. 115551 Martin-Pintado et al. ISUOTj) . 
Such intermittent shocks, created by frequent cloud-cloud 
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collisions can also explain the temperature gradients seen 
in the NH3 data. Cosmic rays can also heat large volumes 
of gas to high temperatures 

It is interesting that the warm gas fraction seems to 
be highest in M 82 where no warm ammonia component 
has been detected. Here the relative ammonia abundance 
in the hot gas component must be < 1CP 9 . This indicates 
that NH3 traces a different molecular gas component in 
M 82 than in NGC 253, Maffei 2 and IC 342. For M 82, the 
presence of an intense dissociating radiation field and the 
predominance of warm molecular filaments with low H2 
column density and density (Mao et al. I2U(JU|) do not only 
explain the low abundance of NH3 (see e.g. the model in 
Fuente et al. I1993[l and other complex molecules but also 
the excitation of CO rotational levels fWei6 l2"001b|) in the 
warm H 2 gas detected by Rigopoulou et al. I|2()()2|l . One 
could argue that larger quantities of warm NH3 exist in 
this low density gas, but are escaping from being detected 
due to lacking excitation. This is, however, improbable 
since NH3 molecules with their low dissociation energy 
would be rapidly destroyed by the interstellar radiation 
field in a low H2 density and column density environment. 

The low NH3 abundance contrasts with a CS abun- 
dance which is normal or even slightly higher than in other 
galaxies. Chemical models of photon dominated regions 
(PDR) show that at moderate extinctions of 6 — 8 mags., 
CS is strongly enhanced up to values of more than 10 -6 
due to the presence of a small amount of ionized sulphur 
(Sternberg & Dalgarno ll995f) . The bulk of the CS we are 
seeing in M 82 might come from such moderately shielded 
regions. 

The lack of ammonia in the warm phase of M 82 also 
suggests that the overall heating mechanism for the molec- 
ular gas is different. In a starburst, photo-dissociating ra- 
diation is capable of heating a large fraction of the gas 
to 150 K (e.g. Rodriguez-Fernandez et al. I2UU1|I . The low 
abundance of ammonia toward M 82, the properties of the 
CO emission (Mao et al. 2000) and the detection of abun- 
dant emission of HCO (Garcia-Burillo et al. l2002f) . a tracer 
of photo-dissociation regions, indicates that PDRs are the 
main heating source that also have influenced the chem- 
istry of the bulk of the circumnuclear gas in this galaxy. 
In regions where the CS abundance is enhanced, shielding 
may be sufficient to prevent the heating of the CS emitting 
gas to high kinetic temperatures. 

5. Conclusions 

The main conclusions of this paper are: 

1. Our detections of highly excited lines of ammonia 
toward the nuclear regions of three nearby galaxies reveal 
the presence of gas components as warm as 50 K . . . 440 K 
toward Maffei 2, IC 342 and NGC 253. Part of the gas may 
be cooler (~ 50 K). These temperatures resemble values 
observed toward the central region of our own Galaxy. 

2. The global beam averaged relative ammonia abun- 
dances toward NGC 253, Maffei 2 and IC 342 are with (3.3 
. . . 5.7) 10~ 8 remarkably similar. The NH3 abundances in 



the warm (<; 150 K) phase in those galaxies, and also in 
the central region of the Milky Way, are ~ 10" 7 . Toward 
M82 this abundance is lower by more than an order of 
magnitude. 

3. Toward NGC 253, IC342 and M82 an elevated H 2 
density of ~ 10 4 cm' 3 is estimated from multilevel studies 
of CS. 

4. The rich chemistry and high ammonia abundances 
observed in NGC 253, Maffei 2, IC342 and the Galactic 
central region may be in part explained by evaporation 
from the ice mantles of interstellar dust grains due to C- 
shocks. These shocks would also provide the heating to 
the observed elevated temperatures. 

5. In the central region of M82, photodissociating ra- 
diation can explain the high temperature of H2 , the lack 
of ammonia and other complex molecules in M 82, and the 
high abundance of CS and HCO. 
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